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Abstract. We present a study of the structural properties of (x)Na2S-(l-x)GeS2 glasses through DFT- 
based molecular dynamics simulations, at different sodium concentrations (0 < x < 0.5). We computed 
the radial pair correlation functions as well as the total and partial structure factors. We also analyzed 
the evolution of the corner- and edge-sharing intertetrahedral links with the sodium concentration and 
show that the sodium ions exclusively destroy the former. With the increase of the sodium concentration 
the "standard" FSDP disappears and a new pre-peak appears in the structure factor which can be traced 
back in the Na-Na partial structure factor. This self organization of the sodium ions is coherent with 
Na-rich zones that we find at high modifier concentration. 



PACS numbers: PACS numbers: 61.43.Bn,6f ,43.Fs,71.15.Pd,63.50.+x 

1. Introduction 

Amorphous (x)Na2S-(l-x)GeS2 systems are interesting glasses, presenting a high ionic conductivity at 
room temperature. These materials can be used as model systems for efficient solid electrolytes even 
though practical applications are more often based on glasses containing lithium or silver ions for 
example. Although the ionic conduction process has been clearly observed in sodium thiogermanate 
glasses, its microscopic origin is still not clearly understood. Like in many other glassy systems, the 
atomic mechanisms responsible for these transport properties have not yet been clearly determined which 
explains the large number of studies dedicated to this topic in the past few years pQ. In particular, it 
is interesting to determine if the modifier ions follow preferential pathways inside the amorphous system 
and if this process is dependent on the glass type. To that end, theoretical studies such as molecular 
dynamics (MD) simulations can provide detailed information concerning the microscopic properties of 
these chalcogenide glasses. 

In previous studies we have analyzed the physical and chemical properties of GeS2 glasses via Density 
Functionnal Theory (DFT)-based molecular dynamics simulations EJ El • The results obtained from 
our model were in very good agreement with the available experimental data. Subsequently we studied 
the vibrational and electronic properties of sodium thiogermanate (Na-Ge-S) glasses using the same 
description 6 . In the present work, we focus on the structural properties of (x)Na2S-(l-x)GeS2 
amorphous systems using the same DFT-based model in order to determine the impact of the Na + ions 
on the structure of germanium disulfide glasses. It should be noted that there is a lack of experimental 
results in the literature concerning these doped chalcogenide glasses. Thus no extensive comparison with 
experiments is possible at this time, and our results should be considered as a first step in the study of 
the structural properties of (x)Na2S-(l-x)GeS2 glasses. 

This paper is organized as follows: section II is dedicated to the description of the theoretical foundations 
of our model. Subsequently, we present the results obtained in our simulations in section III. Finally, 
section IV contains the major conclusions of our work. 



2. Model 



The mathematical model we have used was developed by Sankey and Niklewski 0. It is included in a 
first-principles type molecular-dynamics code, called fireball96, that is based on the DFT within 
the Local Density Approximation (LDA) A tight-binding- like linear combination of pseudoatomic 

orbitals, satisfying the atomic self-consistent Hohcnbcrg-Kohn-Sham equations is used to determine 
the electronic eigenstates of the system. A basis set of one s and three p pseudo-orbitals, slightly excited 
since they vanish outside a cut-off radius of 5a (2.645 A), is required. 



Structural properties of various sodium thiogermanate glasses through DFT-based molecular dynamics simulations'! 



The pseudopotential approximation is used to replace the core electrons with an effective potential which 
acts on the valence electrons, and Hamman-Schluter-Chiang pseudopotentials were used in this work 
|12| . The Hamiltonian is calculated through the use of the Harris functional which avoids long 
self-consistent calculations by a zero*' 1 order approximation of the self-consistent density. Further details 
on the method can be found in the original paper |14| . Periodic boundary conditions are applied in order 
to limit surface effects, and only the T point is used to sample the Brillouin zone. This model has given 
excellent results in several different chalcogenide systems over the past few years El QUI ■ 
In the present work we computationally melted a crystalline a-GeS2 configuration containing 258 particles 
at 2000K during 60 ps (24000 timesteps) until we obtained an equilibrated liquid. Subsequently we 
randomly replaced GeS4 tetrahedral units by Na2S3 units, in order to obtain a given concentration of 
modifier without affecting the global electrical neutrality of the system. This procedure has been used in 
a similar manner in SiC>2 glasses |17|. and keeps the total number of atoms, N, unchanged (258 particles 
in this work). 

Following this scheme, we generated eight (x)Na2S-(l-x)GeS2 samples with different sodium 
concentrations (x= 0, 0.015, 0.03, 0.06, 0.11, 0.2, 0.33 and 0.5). The bounding box was rescaled for each 
sample so that the density matched its experimental counterpart (with cell lengths ranging from 19.21 
A for x=0[IB] to 18.3 A for x=0.5[]2]). Since the experimental densities are not exactly the equilibrium 
densities of the model, the pressure inside the samples is not exactly equal to zero. Nevertheless it 
is extremely small (< 0.4 GPa) demonstrating the quality of the model. The resulting system was 
subsequently melted at 2000K during 60 ps in order for the system to lose memory of the previous 
artificial configuration. When the liquid reached an equilibrium state, the systems were quenched at a 
quenching rate of 6.8x 10 14 K/s, decreasing the temperature to 300K, passing through the glass transition 
temperature T g . The samples were then relaxed at 300K over a period of 100 ps. In order to improve 
the statistics of our results, five samples, starting from independent liquid configurations, were generated 
for each concentration. The results presented here have been averaged over these five samples. 

3. Results 

In a previous work jS| , we determined that the atomic charges in GeS2 glasses using a Lowdin description 
|2Tj] are equal to +0.94 for Ge atoms and -0.47 for the S particles. Even though homopolar bonds exist in 
this glass 0], no Ge-Na bonds were seen in our samples at the concentrations studied here (0 < x < 0.5). 
Therefore we can study the impact of the sodium ions on the short-range order of germanium disulfide 
glasses by analysis of the S-Na bonds exclusively. 
Among the radial pair correlation functions g a f}{f) defined as: 

9a ^ = A^Kdr (1) 
for the different (a,/3) pairs, we can hence concentrate on the gsNair) function. It can be seen in Fig.l that 
the gsNa( r ) for S-Na pairs does not change significantly with the concentration of the sodium modifier. 
The first-neighbor peak appears at 2.4 A for all concentrations (0 < x < 0.5), which is slightly inferior 
to experimental data concerning crystalline sodium thiogermanate systems (2.71 A |21j). However we 
tested several glassy samples with the program siesta |22| in a self-consistent ab initio description, using 
the most accurate basis set available in the code for the description of the atomic orbitals (Double Zeta 
Polarized + Generalized Gradient Approximation): the results were very similar (ds-Na = 2.45 A) to 
those obtained with fireball96. 

The different coordination numbers n a p can also be calculated. They correspond to the average number 
of (3 neighbors of a given a atom within a sphere of radius r m j n , r m j„ being the first minimum of the 
corresponding <? a /3(r). We show in Fig. 2 the evolution of the coordination number of the sulfur atoms 
with the three elements present in our system. It can be seen that the value of ns-Ge decreases with 
the addition of sodium modifier from 1.95 for x=0 to 1.3 for x=0.5. This indicates that the Na + ions 
break several inter-tetrahedral connections in order to be connected to the sulfur atoms, which explains 
the increasing value of ns-Na (2.45 for x=0.5). The value of ns-s relative to homopolar bonds remains 
constant without any dependency with the sodium concentration. Finally the average coordination 
number of the sodium atoms relative to the sulfur atoms, njv a _s, has been found close to 3.27 for x=0.5. 
In fact, the Na + particles in the system are connected to 3 or 4 sulfur atoms, which is consistent with 
the results obtained by Foix et al. in cluster simulations of Na2GeS3 glasses using a Hartree-Fock model. 

eh 

We reported in a previous study £j] that 14.53% of the sulfur atoms in glassy GeS2 systems are non- 
bridging (i.e. connected to only one Ge atom). These bond defects have been found [H] to create a 
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negatively charged electronic environment around the terminal sulfur. Therefore it is interesting to study 
the evolution of these non-bridging sulfur particles with the introduction of Na + cations. We distinguish 
here between the terminal sulfur atoms connected to (a) zero, (b) one, (c) two or (d) three sodium ions, 
and represent in Fig. 3 the evolution of these structural entities in (x)Na2S-(l-x)GeS2 for < x < 0.5. 
One can see that the proportion of terminal sulfur atoms that are not linked to a sodium ion decreases 
significantly as soon as x is non zero. This variation is counterbalanced by the increase of the proportion 
of terminal sulfur atoms connected to one sodium atom, which reaches approximately the same value of 
«15% of the S particles in the sample. We can therefore deduce that, as expected, the non-bridging sulfurs 
present in <?-GeS2 significantly attract the Na + ions. Between x — 0.11 and x — 0.33 these structural 
entities decrease (Fig. 3(b)), and are counterbalanced by the increase of terminal sulfur connected to two 
Na ions (Fig. 3(c)), that also reaches the value of «15 %. Finally, for x > 0.33 the proportion of terminal 
sulfurs connected to three sodium ions reaches «15 % as in the aforementioned cases. This shows that 
the non-bridging sulfur atoms observed in g-GeS2 attract several Na + ions in sodium thiogermanate 
glasses. 

The study of the short-range charge deviation in g-GeS2 has also revealed the existence of "positively" 
charged zones, which are mainly caused by 3-fold coordinated sulfur and homopolar bonds 0. We find 
that these positively charged zones remain basically unchanged in (x)Na2S-(l-x)GeS2 glasses, and that 
the Na + ions never connect to 3-fold coordinated sulfur atoms. Hence negatively and positively charged 
zones in g-GeS2 appear to attract and repel the Na + cations respectively. This observation, which seems 
to be self-evident, is in fact relevant since all the different samples have been melted, quenched, and 
relaxed independently. Despite that, we find similarly charged zones in comparable proportions in all of 
these samples. 

The decrease of the ns-Ge coordination number with increasing sodium content observed in Fig. 2 shows 
that an increasing number of intertetrahedral bonds are destroyed by the addition of sodium ions. Since 
two different types of intertetrahedral links exist in glassy GeS2 (edge-sharing and corner-sharing links), 
it is interesting to determine which type of S-Ge connections are destroyed by the introduction of the Na + 
ions. These two types of intertetrahedral connections can be distinguished in the radial pair correlation 
function of germanium-germanium pairs 0j, leading to Ge-Ge distances equal to 2.91 A (edge-sharing) 
and 3.41 A (corner-sharing). In Fig. 4 we present the total number of corner-sharing and edge-sharing 
links in the different samples as a function of sodium concentration. It can clearly be seen that exclusively 
corner-sharing links are destroyed by the introduction of sodium cations whereas the number of edge- 
sharing links remains constant for all sodium concentrations. This result in not predictable a priori since 
very different results may be obtained depending on the glass and the modifier. For example, experimental 
studies on chalcogenide SiS2 glasses have shown that Li + cations destroy exclusively edge-sharing links 
|23II24) . while the addition of Na + ions was found to break corner-sharing and edge-sharing links in equal 
proportions Here we find that sodium ions destroy exclusively corner-sharing links in GeS2 glasses 
but this result still has to be confirmed (or refuted) experimentally. 

In order to analyze the global impact of the sodium ions on the structure of thiogermanate glasses, we 
studied next the total static structure factor S(q) of our (x)Na2S-(l-x)GeS2 glasses. The total static 
structure factor was computed using the following formulae: 



In the above equations, r max — L/2 where L is the edge of the box, the g a f3{r) are the radial pair 
correlation functions (a, (3 = Ge,S,Na), c a and eg are the concentrations of the species a and /? and b a 
and hp are their scattering lengths taken equal to 8.18 fm for Ge, 2.84 fm for S and to 3.63 fm for Na 26 . 
The results for S(q), calculated using the above definition, are shown in Fig. 5 for the different values of 
x (even though these results are directly comparable to neutron scattering experiments, no experimental 
data is available to our knowledge). 

The first sharp diffraction peak (FSDP), which is a signature of amorphous materials, appears at q =1 
A -1 in <?-GeS2 [3J. This peak reveals the existence of a structural order on a length scale d corresponding 




(2) 



where 




(3) 
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to d = 2n/q rj 6.3 A in germanium disulfide glasses. In sodium thiogermanate glasses, it can be 
seen that the intensity of the FSDP decreases significantly with x for x < 0.2. This decrease of the 
FSDP has also been observed experimentally in silver thiogermanate glasses [21] ■ This means that the 
introduction of ions disturbs the intermediate range order of the glassy matrix by breaking corner-sharing 
intertetrahedral connections. However for x > 0.2, a new feature appears in the structure factor at a 
lower value of q (q = 0.5 A -1 ). For x=0.5 this pre-peak reaches almost the same intensity than the 
original FSDP, which almost completely vanishes and shifts to a slightly higher q value. In order to 
find the origin of this new peak we calculated the partial structure factors (computed from the radial 
pair correlation functions g a p(r)) in the different samples, and in particular the Na-Na partial structure 
factor (Fig. 6). We find that the sodium-sodium pairs are mainly responsible for the pre-peak observed 
in Fig. 5. It can indeed be seen that a peak emerges at q = 0.5 A -1 for x > 0.11 and subsequently grows 
with increasing x. To analyze the structural arrangement of the sodium atoms, which are obviously 
not distributed homogeneously inside the simulation box, we present the density of sodium atoms in a 
(0.33)Na2S-(0.66)GeS2 averaged over the relaxation time in Fig. 7. To that end, we divide the bounding 
box into 1000 (10x10x10) elementary volume elements. We subsequently calculated the number of Na 
atoms located in each specific volume element, and integrated this quantity over the relaxation time. 
Hence each sphere represented in Fig. 7 is located at the center of a volume element, and its size increases 
(becomes more red in color) with the local Na density (the largest sphere and thus the highest density 
corresponds to 0.154 A~ 3 ). This figure shows that our sodium thiogermanate sample contains several 
zones which have a high concentration of Na + and other zones which have a poor concentration of Na + 
ions. When we analyzed the local structure around the sodium atoms we found that most of them are in 
a Na2S-like environment, since for z=0.5, 23.5 % of all the sulfur atoms are exclusively connected to Na 
ions, i.e. without any germanium atom as their nearest-neighbors. These sodium rich zones could be the 
static "remainders" of conduction channels that could be observed at higher temperatures similarly to 
what has been reported in sodium enriched silicate glasses via classical molecular dynamics simulations 
and inelastic neutron scattering experiments [171 127j . If we accept this hypothesis, then, since the new 
pre-peak appears only for high sodium concentrations (see Fig. 6), conduction channels would be absent at 
low sodium content and this could explain the experimentally observed increase in the sodium diffusion 
constant when the sodium concentration is greater than x=0.10 |2Hj- Unfortunately, the DFT-based 
model that we used to perform these simulations are too computationally expensive making it difficult to 
reach the simulation time required to study the diffusion of the Na + ions, and to analyze the evolution 
of these zones from a dynamical point of view. This is the reason why the sodium rich zones presented 
in Fig. 7 cannot be directly assimilated to the diffusion channels. 

Most of theses results have no experimental counterparts and therefore need still to be confirmed or 
invalidated, but neutron diffraction studies on these systems have been performed and the results should 
be available in a not so distant future. 

To pursue this study further it would also be interesting to analyze the influence of the nature of the 
modifier on the results. Thus similar studies on silver thiogermanate glasses are currently in progress. 

4. Conclusion 

We have studied the structural properties of sodium thiogermanate glasses via DFT-based molecular 
dynamics simulations for eight different modifier concentrations. We find that the structure of these 
chalcogenide glasses is significantly influenced by the extended positively and negatively charged zones 
observed in amorphous GeS2 5 . As expected, the sodium cations are mainly attracted by the negative 
zones and avoid the positive regions. At high sodium concentrations, the Na + ions therefore become 
localized in space, and zones containing high and low density of sodium cations can be observed. These 
Na2S-rich zones could be the low temperature analogs of the "conduction channels" observed in silicate 
glasses, where the diffusing ions follow preferential pathways. We have also determined that the addition 
of sodium cations inside the (x)Na2S-(l-x)GeS2 system destroys exclusively the corner-sharing links for 
x < 0.5. 

The structure factor of the different samples shows that the FSDP that appears in amorphous GeS2 
at q =1 A -1 , decreases significantly with the addition of sodium ions. At high sodium concentrations 
(x > 0.11), a new pre-peak appears at q = 0.5A -1 , which corresponds to a distance in real space of 
12.6 A. The Na-Na partial S(q) shows that the Na-Na pairs are responsible of this pre-peak indicating a 
non homogeneous distribution of the Na + ions inside the glass which is consistent with the existence of 
the aforementioned sodium rich zones. 
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Figure 1. Radial S-Na pair distribution functions of (x)Na2S-(l-x)GeS2 glasses for < x < 0.5 




Figure 2. Coordination number of sulfur atoms in sodium thiogermanate samples, plotted against x, 
the concentration of sodium atoms. 
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Figure 3. Proportion of non-bridging sulfur atoms in (x)Na2S-(l-x)GeS2 glasses. We distinguish the 
S atoms connected to (a) zero, (b) one, (c) two and (c) three Na+ ions 
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Figure 4. Number of edge-sharing and corner-sharing links in (x)Na2S-(l-x)GeS2 glasses, plotted 
against x, the concentration of sodium atoms. 
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Figure 5. The simulated static structure factors of the different sodium thiogermanate glasses. 




Figure 6. Partial S(q) for Na-Na pairs at different sodium concentrations versus q. 
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Figure 7. Illustration of the local Na density in a (0.33)Na2-(0.66)GeS2 sample integrated over the 
relaxation time (100 ps). The largest sphere corresponds to a Na density of 0.154 A -3 . 



